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Abstract 

The spatiotemporal nonlinear Schrodinger equation with power-law nonlinearity in PT-sym metric potentials is 
investigated, and two families of analytical three-dimensional spatiotemporal structure solutions are obtained. The stability 
of these solutions is tested by the linear stability analysis and the direct numerical simulation. Results indicate that solutions 
are stable below some thresholds for the imaginary part of PT -symmetric potentials in the self-focusing medium, while 
they are always unstable for all parameters in the self-defocusing medium. Moreover, some dynamical properties of these 
solutions are discussed, such as the phase switch, power and transverse power-flow density. The span of phase switch 
gradually enlarges with the decrease of the competing parameter k in PT -symmetric potentials. The power and power-flow 
density are all positive, which implies that the power flow and exchange from the gain toward the loss domains in the VT 
cell. 
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Introduction 

In the last few decades, there has been a surge of interest in 
obtaining exact analytical solutions of nonlinear partial differential 
equations (NPDEs) to describe the natural physical phenomena in 
numerous branches from mathematical physics, engineering scienc- 
es, chemistry to biology [1-3]. Exact solutions often facilitate the 
testing of numerical solvers as well as aiding in the stability analysis. 

The nonlinear Schrodinger equation (NLSE), as one of important 
nonlinear models, has now become an intensely studied subjects due 
to its potential applications in physics, biology and other fields. 
Abundant mathematical solutions and physical localized structures 
for various NLSEs have been reported. For example, bright and 
dark solitons and similaritons [4-6], rogue waves [7], nonautono- 
mous solitons [8] and light bullets [9] etc. have been predicted 
theoretically and observed experimentally in different domains. 

Recently, two-dimensional accessible solitons [10] and nonau- 
tonomous solitons [11] for NLSE in parity-time (VT)- symmetric 
potentials have been reported. The PT-symmetry originates from 
quantum mechanics [12], and was introduced into optical field 
since the important development on the application of 
PTsymmetry in optics was initiated by the key contributions of 
Christodoulides and co-workers [13]. Quite recently, various 
nonlinear localized structures in PT-symmetric potentials have 
been extensively studied. Nonlinear localized modes in PT- 
symmetric optical media with competing gain and loss were 
studied [14]. The dynamical behaviors of (l + l)-dimensional 



solitons in PT-symmetric potential with competing nonlinearity 
were investigated [15]. Bright spatial solitons in Kerr media with 
PT-symmetric potentials have also been reported [16]. Dark 
solitons and vortices in PT-symmetric nonlinear media were 
discussed, too [17]. Moreover, Ruter et al. [18] and Guo et al. [19] 
studied the experimental realizations of such VT systems. 
However, three-dimensional (3D) spatiotemporal structures in 
PT-symmetric potentials are less studied. Especially, 3D spatio- 
temporal structures in PT-symmetric potentials with power-law 
nonlinearities are hardly reported. 

The aim of this paper is to present 3D spatiotemporal structures 
of 3DNLSE with power-law nonlinearity in PT-symmetric 
potentials. Two issues are firstly investigated in this present paper: 
i) analytical spatiotemporal structure solutions are firstly reported 
in PT-symmetric power-law nonlinear media, and ii) linear 
stability analysis for exact solutions and direct simulation are firstly 
carried out in PT-symmetric power-law nonlinear media. Our 
results will rich the localized structures of NLSE in the field of 
mathematical physics, and might also provide useful information 
for potential applications of synthetic PT-symmetric systems in 
nonlinear optics and condensed matter physics. 

Results 

Analytical spatiotemporal structure solutions 

The propagation of spatiotemporal structures in a PT- 
symmetric nonlinear medium of non-Kerr index is governed by 
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the following 3DNLSE 

iu z + PiA L u + p 2 u tt + yJu\ 2m u+[V(x)ViW(x)]u = ^, (1) 

where Aj_=(d 2 ,d 2 ),x = (ri,r 2 ,r 3 ) = (x,y,t), the complex envelope 

of the electrical field u(z,r) is normalized as (koWo)~ l (n 2 /no)~ 1 ^ 2 
with linear index n 0 and Kerr index n 2 , longitudinal z, transverse 
coordinates and comoving time t are respectively scaled to the 
diffraction length hjy=k^w\, the input width unit wq with the 
wavenumber k$ = Ikhq / X at the input wavelength X and \fLj). 
Parameters Pi and P 2 are respectively the coefficients of the 
diffraction and dispersion, and y m for m=l,...,n stand for the 
nonlinearities of orders up to 2/z+ 1. For m—\ one has the 
simple Kerr nonlinearity, for m — 2 the quintic, for m — 3 the septic, 
and so on. Functions V(x) = k^wldnR(x) and W(x) = k^wldni(x), 
with the perturbation of index by a complex profile 
n = no[\ +«/{(r) + iw/(r)], are the real and imaginary components 
of the complex 'PT-symmetric potential, and correspond to the 
index guiding and the gain or loss distribution of the optical 
potential respectively. V and W satisfy V(x) = V( — r) and 
W(x)=-W(-x). 
We consider that solutions of 3DNLSE (1) is of the form: 



u(x) = <I>(r) exp [i/iz + i® (r)] , 



(2) 



where two real valued functions and 9 satisfy the following 
differential equations: 



j2Bj(Aj^-\Vj@\ 2 (S>)+[V(r)-n](S>+y m ^ 2m+1 = 0 (3) 

7 = 1 



J2 Bj(<S>Aj® + 2Vj©Vj<S>) + W(r)Q> = 0. 



(4) 



where = sinh(r y ) and 2 Fi(a,b,c,o) is the Hypergeometric 
function [20]. The parameters in the potential (4) and (5) satisfy 
that V n = K 12 =2ft,K 13 =2feF 21 = V 22 = W 2 /[(k + 2) 2 p x ], V 23 
= W 2 /[(k + 2) 2 p 2 \,W i = W 2 ,n = 2p i +P 2 , with three arbitrary 
constants V 3 , W\ and W$. Parameter m has a serious impact on 
the nature of the gain and loss profile W{r). The value of k as zero 
or nonzero leads to Wir) as asymptotically non-vanishing or 
localized (asymptotically vanishing), respectively. For instance, if 
k = 0, it is the first type of extended Rosen-Morse potential, and if 
k = 1 , it is the first type of extended hyperbolic Scarf potential. 

From (6), V 3 y m <0, thus solution (6) exists in self-focusing (SF) 
media with positive nonlinearity (y m > 0) if V 3 < 0, as well as in self- 
defocusing (SD) media with negative nonlinearity (y m < 0) if 
V 3 >0. 

Specially, if the value of A: is chosen as 0-3, 0(r) has the different 
forms shown in Table 1. 

Case 2 Second type of extended PT-symmetric 
potential. In the following PT-symmetric potential 



V{x)= [^iysech 2 (r 7 -)+ V 2j sech 2k (rj)} + V 3 ft sech 2 (r y ), 



7 = 1 
3 



(7) 



W(x) = ^sechV 7 -)tanh(r/), 

7=1 



with real parameters Vy, Vy, V 3 and Wj, and the competing 
parameter k, the localization condition 0^0 as r-^±oo leads to 
solution of Eqs. (3) and (4) in the form 



0>(r) = 



V 3 



sech(x)sech(y)sech(r) 



(8) 



7=1 



mWj 

H(r)= 4^(km + 2)Bj 



2F1 



1 £+1 3 



in 



where j = 1 ,2,3 A = B 2 = ft ,£ 3 = ftA = d 2 , A 2 = d 2 ,A 3 = d 2 t ,V x 

= d x ,v 2 = d y y 3 = d t . 

In the following, we obtain analytical spatiotemporal structure 
solutions of Eqs. (3) and (4) in two different PT-symmetric 
potentials. 

Case 1 First type of extended PT-symmetric 

potential. Considering the PT-symmetric potential 



3 

V(x) = Y [Kysech 2 (#)-) + V 2j SQch 2k (r y )] + V 3 Ii sech 2m (r y ), 
j=i " ' J = i 

3 

W(x) = Y F^sech^)tanh(r y -), 

7 = 1 



(5) 



with real parameters Vy,V 2 j,V 3 and Wj, and the competing 
parameter k, the localization condition O— >0 as r^±oo yields 
solution of Eqs. (3) and (4) 



(D(r) = 



V 3 



2m 



sech(x) sech(y) sech(0 , 



®(r)=E 



Wi 



1 *+l 3 _ 2 

2' 2 '2' y 



(6) 



with constant V\\ = V\ 2 = (m+\)P X / m 2 , V\ 3 = (m+\)P 2 / m 2 , 
V 2l = V 22 = m 2 W 2 [/ [(km + 2) 2 ft ] , K 23 = m 2 P^ 2 / [(km + 2) 2 ft ] , 
J^2 = W\,n = (2 fa + ft)/m 2 and three arbitrary constants V 3 ,Wi, 
and W^. Moreover, <^ = sinh(r y ) and 2 F\(a,b,c,o) is the Hyper- 
geometric function. 

Specially, when the value of k is chosen as 0-3, the expressions 
of 9(r) are shown in Table 2. 

Note that for the Kerr nonlinearity (i.e. m = 1), solutions (6) and 
(8) are same. 

Properties of spatiotemporal structure solutions 

The even and odd functions for the real part V and imaginary 
part W of the PT-symmetric potential (7) are shown in Fig. 1 in 
regard to x^y and t for different L Figs. 1(c) and (d) show V for 
different k at z = 30, y = 0, t=\0 when m — 2 and 1, respectively. 
Fig. 1(e) shows Wfor different k at z = 30, y = 0, t = 10 when m = 2. 
From the yellow dash lines in Figs. 1(c)— (e), V is localized when 
m = 1 or 2, while W is asymptotically non-vanishing in the 2D 
extended Rosen-Morse potential. It possesses unbroken VT- 
symmetry [18]. From red crosses, blue lines and black circles in 
Figs. 1(c)— (e), the peaks and widths of Fand ^gradually decrease 
when k increases. Compared red crosses, blue lines and black 
circles in Fig. 1(c) with those in Fig. 1(d), the amplitudes of V 
attenuate when m adds. 
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Table 1. The expression of ®(r) in first type of extended ;PT-sym metric potential. 



0(r) 

— — { arctan [sinh(.x)] + arctan [sinh(y)j } + — — arctan [sinh(7)] 
3pi 3p 2 

W\ Wi 

— - [tanh(x) + tanh(y)j + — -tanh(7) 



Wi 



{ arctan [srnh(x)] + tanh(x)sech(x) + arctan [sinh(y)] + tanh(y)sech(y)} + A { arctan [sinh(7)] + tanh(7)sech(7)} 



doi:1 0.1 371/journal.pone.01 00484.t001 



In the PT-symmetric potentials above, we can find the 
phenomena of phase switch of solutions (6) and (8). Fig. 2(a) 
exhibits switches of phase in (8) for different PT-symmetric 
potentials (7) with m = 2. When k decreases, phases switch from 
smaller to bigger values along x, and the spans of switch gradually 
enlarge. However, in the Rosen-Morse potential with k=0, no 
phase switch appears. In the 'PT-symmetric potential (5) with 
m — 1,2, the phase switch also exists. We omit the related plots. 

The power P can be expressed as P = |w(z,r)| 2 
dr = 8(— Vi/y m ) l / m for solution (6). For solution (8), when m = 1, 
P=-8V 3 /y u and when m=2, P= -n 3 (- V 3 /y 2 ) l/2 . These 
powers are all positive for any values of the parameter V$ due to 
V 3 y m <0, and independent of the parameter of the imaginary part 
in the PT-symmetric potential. Moreover, the power-flow 
densities (Poynting vectors) S = (i/2)[OVO* — €>* V€>] across spa- 
tiotemporal structure solutions (6) and (8) have the form S = 
Is-Kij^ SQC h k+2 (x), f s Q ch k+2 (y), f sech* +2 (0] 



s 



and 



(- n)*[m sec h*+*(x),^sech* + %), 



f sech^O], 



W\ 

which are related to the competing parameter k and parameter 
m. Obviously, due to V 3 y m <0, S is everywhere positive, which 
indicates that the power flow and exchange for solutions (6) and (8) 
in the VT cells are always from the gain toward the loss domains 
(one direction). An example to this case is shown in Fig. 2(b) for 
k=3. The similar results also exist when other k and m are chosen. 

Discussion and Analysis 

Linear stability analysis of analytical solutions 

We study the linear stability of solutions (2) with (6) and (8) 
via the method developed in [21] when m— 1,2. A perturbation 
of an exact solution can be expressed as 



u(z,r) = {w„(r) + £[i?(r) + /(r)]exp(iaz)}exp(i^z), where e is an 
infinitesimal amplitude, u n (r) is a solution of Eq.(l), R(r) and I(r) 
are the real and imaginary parts of perturbation solution, and o 
represents the perturbation growth rate. Inserting this expression 
into Eq. (1) and linearizing it around the unperturbed one (the 
first-order term of e), we have the eigenvalue problem 



L + R = al, 
L_I = oR, 



(9) 



where o is an eigenvalue, R and / are eigenfunctions with 
Hermitian operators L+ = — P\{d 2 x + dy) — fad 1 — f] + yiU n (r) 2 
— {V + iW) + n with rj+ = 3 and rj- = 1 for m = 1 and 
L ± = -fa{d 2 x + d 2 y )-fad 2 t -v ± y 2 u n {x) 2 -{V + iW) + ii with v + 
= 5 and v_ = 1 for m = 2. If all imaginary parts of a are equal to 
zero, solution can be completely stable. Otherwise, if any 
eigenvalue a possesses an imaginary part, the perturbed solution 
would add exponentially with z and thus corresponding solution 
becomes linearly unstable. 

The eigenvalues a of solutions (6) and (8) in the SF and DF 
media under the 2D extended Rosen-Morse potential have many 
imaginary parts, and thus solutions (6) and (8) are always unstable 
in these nonlinear media. Fig. 3 shows some examples of the 
eigenvalue o in the SF and DF media. From Figs. 3(a) and 3(b), the 
eigenvalues a for both SF and SD nonlinearities have many 
imaginary parts, and thus solutions (6) and (8) with m = 1 are 
unstable. Similar, solutions (6) and (8) with m = 2 are also unstable 
because there exist many imaginary parts of the eigenvalue a in 
Figs. 3(c) and 3(d), too. The asymptotically non-vanishing 
characteristic of W in the 2D extended Rosen-Morse potential 



Table 2. The expression of 0(r) in second type of extended VT -symmetric potential. 



0(r) 



mWx , , mWi 

YYl Wx WI W~2 

{ arctan [sinh(x)] + arctan [sinh(y)j } + — — arctan [sinh(f)] 



(m + 2)£ 
m W\ 
(2m + 2)/? : 

mWx 



(m + 2)P 2 



[tanh(x) + tanh(y)j- 



(2m + 2)# 



■tanh(0 



— — _|_2)ff { arctan [sinh(x)] + tanh(x)sech(x) + arctan [sinh(y)j + tanh(y)sech(y)} + 2(3^+2)^ {arctan [sinh(?)j + tanh(7)sech(/)} 
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-5 -5 





Figure 1. The VI -symmetric potential (7): (a) and (b) Isosurface plots of V and W for k = 

z = 30, y = 0, t= 10 when m=2 and 1, respectively; (e) 1/1/ for different k at z = 30, y = 0, t 

ft = 1.1,0 2 = 1.2, ^i=0.8, ^3=0.9, K 3 = -13. 
doi:1 0.1 371 /journal. pone.0100484.g001 



3 at z = 20; (c) and (d) 1/ for different k at 

= 10 when m = 2. Parameters are chosen as 



leads to the linear instability of solutions (6) and (8) with m = 1 and 
m=2. 

Under the 2D extended hyperbolic Scarf potential, solutions (6) 
and (8) with m = 1 and m = 2 are stable below some thresholds for 
W\ and W 3 in the SF medium, while they are always unstable for 
all parameters in the SD medium. Fig. 4 exhibits some examples of 
the eigenvalues a in the SF and DF media. From Figs. 4(a), (c) and 
(e), the eigenvalues er of solutions (6) and (8) with m — 1 and m = 2 



are all real, and thus solutions are linearly stable in the SF 
medium. When ft = 1.1, fi 2 = 1.2, V 3 = - 13, k= 1, y x = 1.1 or 
y 2 = 1-1, the thresholds are W\ <0.04, W3 < 0.043 for solutions (6) 
and (8) with m= 1, ^ < 0.051, J^ 3 < 0.056 for solution (6) with 
m = 2, and fl^i < 0.005, ^3 < 0.006 for solution (8) with m=2, 
respectively. However, solutions (6) and (8) with m— 1,2 are always 
unstable in the SD medium because there exist some imaginary 
parts of the eigenvalues o for all parameters. Some cases are shown 




Figure 2. Switches of phase (8) for different VT -symmetric potentials (7) in (a). Power-flow vector S for solution (8) when k = 3 indicating 
the power flow from gain towards loss domains in (b). Parameters are chosen as the same as those in Fig. 1. 
doi:1 0.1 371 /journal. pone.01 00484.g002 
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x 10 



(a) 



¥ 0 



-2 



Case 1,2an8l,kaO,Y <| 8l.1 



:xa aa<t 



x 10 



(b) 



E 0 



Case 1,2:m=1 ,k=0,y^=- 1.1 




-2 0 

Re(a) 

3 (d) 



x 10 



Case 2:m=2,k=0,7 2 =1 -1 



-2 0 2 

Re(o) 

Figure 3. Eigenvalues for solution (6) and (8) in (a),(c),(d) SF medium and (b) SD medium under the 2D extended Rosen-Morse 
potential. Parameters are chosen as W\ =0.04, JV 3 =0.043, ft = 1.1, ft = 1.2 with (a),(c),(d) y 3 = -13 and (b) V 3 ='\3. Other parameters are shown in 
the plots. 

doi:1 0.1 371 /journal.pone.01 00484.g003 



in Figs. 4(b), (d), (f). From these results, the gain (loss) related to the 
values of W\ , W 3 should be enough small compared with a fixed 
value of V 3 , otherwise, solutions (6) and (8) with m = 1 and m = 2 
eventually lead to instability. 

Furthermore, when £=2,3 in the 2D extended 'PT-symmetric 
potentials (5) and (7), solutions (6) and (8) with m = 1 and m = 2 are 
stable below some thresholds for W\ and W 3 in the SF medium 
because the eigenvalues a of solutions (6) and (8) with m = 1 and 
m = 2 are all real from Fig. 3. When ft = 1.1, ft = 1.2, 
F3 = — 13, k = 2, y x = 1.1 or y 2 = l-l 5 tne thresholds are 
W\ < 0.043, FF 3 < 0.051 for solutions (6) and (8) with m= 1, 
W 7 ! < 0.06, < 0.066 for solution (6) with m=2, and 
W\ < 0.0072, P^ 3 < 0.058 for solution (8) with m = 2 from 



Figs. 5(a), (c), (e), respectively. For 



= 1.1, ft = 1.2, V 3 



-13, 



k = 3, 7i = 1.1 or y 2 = l.l, the thresholds are PTi < 0.047, 
P^3< 0.057 for solutions (6) and (8) with m=l, W\< 
0.074, W 3 < 0.082 for solution (6) with m = 2, and ^ = 
W 3 < 0.0061 for solution (8) with m=2 from Figs. 5(b), (d), (f), 
respectively. However, in the SD medium, solutions (6) and (8) 
with m— 1,2 are always unstable because there also exist some 
imaginary parts of the eigenvalues a for all parameters. 

When k is bigger, we have the similar results. Solutions (6) and 
(8) with m — 1 and m = 2 are stable below some thresholds for W\ 
and W$ in the SF medium, while they are always unstable for all 
parameters in the SD medium. Here we omit these discussions. 

Numerical rerun of analytical solutions 

Based on the linear stability analysis, we know the stable 
domains of analytical solutions under different 2D extended VT- 
symmetric potentials. In the following, we further test the stability 
of these solutions by the direct numerical simulation. Here we use 
a split-step Fourier pulse technique. In real application, the 
analytical cases are not exactly satisfied, thus we consider the 
stability of solutions with respect to finite perturbations. The 
perturbations of 5% white noise are added to initial fields coming 
from solutions (6) and (8) of Eq. (1). 
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.1.1 
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2 



-202 
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Figure 4. Eigenvalues for solution (6) and (8) in (a),(c),(e) SF 
medium and (b), (d), (f) SD medium under the 2D extended 
hyperbolic Scarf potential. Parameters are chosen as ft = 1.1, 
ft = 1.2 and (a),(b) Parameters are chosen as ft = 1.1, ft = 1.2 and 
(a),(b) ^1=0.04,^ = 0.043, (c),(d) W x =0.051, W 3 =0.056, (e),(f) 
^1=0.005,^3=0.006 with (a),(c),(e) l/ 3 =-13 and (b), (d), (f) V 3 = -\3. 
Other parameters are shown in the plots. 
doi:1 0.1 371/journal.pone.01 00484.g004 



Figure 6 exhibits the numerical reruns corresponding to 
Figs. 4(a)-(f) in the 2D extended hyperbolic Scarf potential. In 
the SF medium, the single VT complex potential is strong enough 
to suppress the collapse of localized solutions caused by diffraction, 
dispersion and different nonlinearities. The numerical solutions in 
Figs. 6(b), (d) and (f) do not yield any visible instability, and good 
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-6 (d) 
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Re(a) 

7 (f) 
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-1 



Case 2:m=2,k=3,Y 2 =1.1 



-2 0 2 

Re(a) 

Figure 5. Eigenvalues for solution (6) and (8) in the SF medium 
under the 2D extended VT -symmetric potential. Parameters are 



chosen as ft = 1.1, 



= 1.2, Vx 



-13 with (a) =0.043, 1^ = 0.051, 



(b) ^i=0.054, ^3=0.057, (c) ^i=0.06, ^ 3 = 0.066, (d) W x =0.074, 
^3=0.082, (e) ^1=0.0072,^=0.058 and (f) ^ = ^ = 0.0061. 
Other parameters are shown in the plots. 
doi:10.1 371 /journal. pone.01 00484.g005 

agreement with results from the linear stability analysis for 
analytical solutions is observed. Numerical calculations indicate 
no collapse, and stable propagations over tens of diffraction/ 
dispersion lengths are observed except for some small oscillations. 
Compared Fig. 6(b) with Fig. 6 (d) or Fig. 6 (f) respectively, we see 
that for the same solution (6) or (8), solution with m — 1 is more 
stable than solution with m = 2 because there are smaller 
oscillations in Fig. 6(b) than those in Fig. 6 (d) or Fig. 6 (f). In 
the DF medium, solutions (6) and (8) are both unstable in the 2D 
extended hyperbolic Scarf potential, which is shown in Figs. 6(c), 



(e) and (g). They can not maintain their original shapes, change 
from distortion to collapse, and ultimately decay into noise. 

Figure 7 displays other examples of stable analytical solutions, 
and it is the numerical reruns corresponding to Figs. 5(a),(b),(d),(f) 
in the 2D extended 'PT-symmetric potential. In the SF medium, 
we can obtain stable spatiotemporal structures. From 
Figs. 5(a),(b),(c),(e), the influence of initial 5% white noise is 
suppressed, and these spatiotemporal structures (6) and (8) stably 
propagate over tens of diffraction/dispersion lengths and only 
some small oscillations appear when k is chosen 2 or 3 in the 2D 
extended PT-symmetric potential. However, in the DF medium, 
spatiotemporal structures are unstable and broken down propa- 
gating after tens of diffraction/ dispersion lengths, and at last turn 
into noise. Compared Fig. 6(d) with Fig. 7 (c) or Fig. 6 (f) with Fig. 7 
(e) respectively, spatiotemporal structures are more stable in the 
2D extended T^T-symmetric potential with k — 3 than those with 
k= 1. 

Conclusions 

We conclude the main points offered in this paper: 

• Analytical spatiotemporal structure solutions are firstly report- 
ed in T^T-symmetric power-law nonlinear media. 

We obtain two families of analytical three-dimensional 
spatiotemporal structure solutions of a spatiotemporal NLSE 
with power-law nonlinearity in T^T-symmetric potentials. 
Some dynamical characteristics of these solutions are dis- 
cussed, such as the phase switch, power and power-flow 
density. The spans of phase switch gradually enlarge with the 
decrease of the competing parameter k in T^T-symmetric 
potentials. The power and power-flow density are all positive, 
which implies that the power flow and exchange from the gain 
toward the loss domain in the VTcell. 

• Linear stability analysis for exact solutions and direct 
simulation are firstly carried out in T^T-symmetric power-law 
nonlinear media. 




Figure 6. Initial value of solution (6) at z = 0 in (a); (b)-(g) the numerical reruns corresponding to Figs. 4(a)-(f) in the 2D extended 
hyperbolic Scarf potential at z = 80. An added 5% white noise are added to the initial values. All parameters are chosen as the same as those in 
Fig. 4. 

doi:1 0.1 371 /journal. pone.01 00484.g006 
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Figure 7. The numerical reruns corresponding to Figs. 5(a),(b),(d),(f) in the 2D extended PT-symmetric potential at z = 80 in 
(a),(b),(c),(e). (d) and (f) are corresponding to (c) and (e) in the SD medium with = 13,y 2 = — 1.1. An added 5% white noise are added to the initial 
values. All other parameters are chosen as the same as those in Fig. 5. 
doi:1 0.1 371 /journal.pone.01 00484.g007 



The stability of exact solutions is tested by the linear stability 
analysis and the direct numerical simulation. Results indicate 
that solutions are stable below some thresholds for the 
imaginary part W of PT-symmetric potentials in the SF 
medium, while they are always unstable for all parameters in 
the SD medium. 

Our results will rich the localized structures of NLSE in the 
field of mathematical physics, and might also provide useful 
information for potential applications of synthetic PT- 
symmetric systems in nonlinear optics and condensed matter 
physics. 
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